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Abstract 

The photo-production of J/'^-^He bound state He\ji^) on a ^He target has been inves- 
tigated using the impulse approximation. The calculations have been performed using several 
7 -|- — )• J/^' + models based on the Pomeron-exchange and accounting for the pion-exchange 
mechanism at low energies. The J/^' wavefunctions in pi?e]j/,j, are generated from various J/\I'- 
nucleus potentials which are constructed by either using a procedure based on the Pomeron-quark 
coupling mechanism or folding a J/^-N potential (wj/^_jv) ™to the nuclear densities. We consider 
uj/ij, TV derived from the effective field theory approach, Lattice QCD, and Pomeron-quark coupling 
mechanism. The upper bound of the predicted total cross sections is about 0.1 — 0.3 pico-barn. We 
also consider the possibility of photo-production of a six quark-J/^ bound state ([^^Jj/^i) on the 
^He target. The Compound Bag Model of A^A^ scattering and the quark cluster model of nuclei 
are used to estimate the [(?^]-N wavefunction in ^He by imposing the condition that the calculated 
^He charge form factor must be consistent with what is predicted by the conventional nuclear 
model. The upper bound of the predicted total cross sections of 7 +^ He — )• [(/^] j/iji + A^ is about 2 
- 4 pico-barn, depending on the model of7 + A^— )■ J + N used in the calculations. Our results 
call for the need of precise measurements of7+p— )• J/^'+p and also the 7 +^ — )• J/^' + n + p 
reactions near the threshold. 

PACS numbers: 25.20.Lj, 24.85. +p 
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I. INTRODUCTION 



The role of the gluon field in determining the interactions between nucleons and quark- 
antiquark (QQ) systems, which do not share the same up and down quarks with the nucleon, 
is one of the interesting subjects in understanding Quantum Chromodynamics(QCD). An 
important step toward this direction was taken by Peskin[lt] who applied the methodology 
of the operator product expansion to evaluate the strength of the color field emitted by 
heavy QQ systems. His results suggested [2] that the van der Waals force induced by the 
color field of J/\l/ on nucleons can generate an attractive short-range J/\l/-iV interaction. By 
using the effective field theory method, Luke, Manohar, and Savagejl] used the results from 
Peskin to predict the J/\l/-nucleon forward scattering amplitude which was used to get an 
estimation that J/\l/ can have a few MeV/nucleon attraction in nuclear matter. Brodsky and 
Miller further investigated the J/\l/-N forward scattering amplitude of Ref.j^ to derive a 
J/^-N potential (f j/ij-.a?) which gives an J/\l/-N scattering length of -0.24 fm. The result of 
Peskin was also used by Kaidalov and Volkovitskyjoj], who differed from Ref. [3] in evaluating 
the gluon content in the nucleon, to give a much smaller scattering length of -0.05 fm. In 
a Lattice QCD calculation, Kawanai and Saski[5| obtained an attractive J/\l/-N potential 
i^j/^,N = - <ye~^^/r with a = 0.1 and fi = 0.6 GeV, which gives a scattering length - 
0.09 fm. In Ref. 0], Brodsky, Schmidt and de Teramond proposed an approach to calculate 
the potential between a cc meson and a nucleus by using the Pomeon-exchange model of 
Dannachie and Lanshoff[8|. The J/\l/-N potential obtained in this approach is fj/iii,7v = - 
ae~^'^ /r with a = 0.6 and /i = 0.6 GeV which gives a rather large scattering length - 8.83 
fm. 

Our first objective in this paper is to explore whether these J/\l/-N potentials, with 
rather different attractive strengths, can form J/\l/-nucleus bound states. Following the well 
developed method in nuclear reaction theory jol], this is done by searching for bound states 
by solving the Schrodinger equation with a folding potential constructed by integrating the 
J/\l/-N potential over the nuclear density. We will also consider the approach of Ref. 0] 
in predicting J/\l/-nucleus bound states by the coherent sum of Pomeon-exchange between 
quarks in J/\l/ and all quarks in the nucleus. For each of the predicted bound pifejj/^, 
systems, we then estimate the photo-production cross section of the ■j+^He — )■ fHe]j/<i, + N 



reaction to facilitate future experimental investigations [10 



The second part of this work is motivated by the investigations by Brodsky and de 
Teramond ll| who found that the spin correlation of pp elastic scattering near the J/\& pro- 
duction threshold can be explained if one postulates the excitation of a hidden charm (c) 
state \qqqqqqcc >. Based on the similar consideration on the role of multi-quark configura- 
tions, Brodsky, Chudakov, Hoyer, and Laget[ll] suggested in a study of — )■ J/'^+n+p 
reaction that J/\I^ can interact strongly with the six-quark [q^] component of the deuteron 
wavefunction because the octet 3-quark [g^Jg in the [q^] can directly interact with each quark 
in J/\I^. These two works suggest the possibility that if J/\l/ overlap with a [q^] cluster in 
nuclei, a bound system could be formed. It is of course very difficult, if not impos- 

sible, to estimate [g^]- J/\l/ interaction. Instead, we will simply assume the existence of such 
states and use the previous works on quark clusters in nuclei to explore how the cross 

sections of 7 -|-^ He — )■ [g^lj/^- + N depend on the parameters characterizing the [g^j-J/^' 
interaction within a potential model. 

Our first task is to construct a model of --f + N ^ J /'^ + N reaction. At high energies, it 
is well recognized that this reaction can be described by the Pomeron-exchange model with 
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an interpretation jsl, [21-23| that Pomeron-exchange is due to the exchange of gluons within 
QCD. This is illustrated in part (a) of Fig{T] At low energies, one expects that mechanisms 
other than Pomeron-exchange could also contribute as can be seen in the exclusive photo- 



production reaction on the nucleon|23l-l25 



However, very little investigation has been done 
for J/\E' photo-production in the near threshold region. As a first step, we will only consider 
the meson-exchange mechanism which can be calculated from using the partial decay width 
of J/\i/ — )■ vrp listed by Particle Data Group (26| (PDG). With the vector meson dominance 
(VDM) assumption, this observed decay process indicates that J /'^ photo-production can 
also be due to the exchanges of a vr meson with the nucleon, as illustrated in part (b) of 
Figffl 
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FIG. 1: Reaction mechanisms of — )• J/^ + : (a) Pomeron-exchange, (b) pion-exchange. 

We next consider the photo-nuclear reaction mechanism that a J/^' is produced from a 
nucleon in a nucleus with mass number A and then forms a bound state with the spectator 
B system which can he a {A — 1) nuclear system or a quark cluster [g^*^"^"^-'] in the target 
nucleus A. With this commonly used impulse approximation, the reaction cross sections 
can be calculated from the 7 + )■ J /'^ + N amplitude, which will be generated from the 
Pomeron-exchange and pion-exchange mechanisms described above, and the initial nucleon 
and final J/\l/ wavef unctions. The nucleon wavefunctions can be taken from the available 
nuclear models. The J/\E' wavefunctions will be generated from various J /^!-B potentials 
mentioned above. For simplicity, we only present the predictions of the cross sections of 
7 + [^He] N + [^He\ji^ and 7 + ^He] N + {q^]j/^ reactions. 

In section II, we present formula for calculating the ''y + N^J/'^ + N amplitudes from 
the Pomeron-exchange and pion-exchange mechanisms. The impulse approximation formula 
for calculating the cross sections of 7-I- [A\ — )• [B] j/^ + are given in section III. Our results 
are presented in section IV. In section V, we give a summary and discuss the necessary 
future work. 



II. FORMULA FOR 7 + iV ^ J/^ + N REACTION 

In the center of mass system, the differential cross section of 7(g) -|- N{—q) — )• J /'^{k) + 
N{—k) with invariant mass W can be written 

^ = [ ^ ][ ^ ]l 2^ 2^ \ < k\j/^mMW)\q\^m, > \ , (1) 

where Xj/^, and A-y are the helicities of the J/\i/ and photon, respectively, rris is the z- 
component of the nucleon spin, and Ea{p) = [m^ + p'^Y^'^ is the energy of a particle with 
mas nia- The reaction amplitude is written as 



rx /I /TT^M^A 1 1 / niN I rriN 1 

< kAj/if,m^\t[W)\qXyms > 



(27r)3 y/2Ej/^{k)\l EN{k) V EN{q) 
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x[u^dPy,{k,\j/^)M^^''{v\v)^u{qA,)urnXp)], (2) 

where Umsip) is the nucleoli spinor (with the normahztion Um,{p)um.'^(j)) = Sms,m'J , 
efj,{k, Aj/*) and e^(g, A^) are the polarization vectors of J/\l/ and photon, respectively. Here 
we also have introduced the four-momenta for the initial and final nucleons: 

p = {EN{q),-q) ■ p' = iE^ik),-k). 

In the following subsections, we give formula for calculating the invariant amplitude M'^'' 
due to the Pomeron-exchange and meson-exchange mechanisms, as illustrated in FiglH 



A. Pomeron-exchange amplitude 



Within the Pomeron-exchange model of Donnachie and Landshoff [8|], the vector meson 
photo-production at high energies is due to the mechanism that the incoming photon couples 
with a qq pair which interacts with the nucleon by the Pomeron exchange before forming 
the outgoing vector meson. The quark-Pomeron vertex is obtained by the Pomeron-photon 
analogy [si, which treats the Pomeron as a C = +1 isoscalar photon, as suggested by a 

Following the formula given explicitly 

(3) 



study of non perturbative two-gluon exchanges [21 
in Ref . [13] , we then have 

M^/{p',p)=GAs,t)rr{p',p) 
with 



mp',p) = ^l'2^/^e 



fy M^-t\2fi', + M^-t 



F,{tmg^^-k^Y] 



(4) 



where t = {p — p'Y, s = {q + pY = W^, = e^/47r, f3g is the Pomeron-quark coupling 
constant. My is the vector meson mass, and Fi{t) is the isoscalar electromagnetic form 
factor of the nucleon, 



AMI - 2-8t 



(4M2,-t)(l-t/0.71)2- 

Here t is in unit of GeV^, and Mn is the proton mass. 
The Regge propagator for the Pomeron in Eq. ([3]) is 



(5) 



Gm — I 



exp 



where a 
In Eq 



pit) 



ao + Ci'pt- It is common j27l| to use ao = 1.08 and a'p 
fy is the vector meson decay constant: = 5.33, - 



= 1/^0 
15.2, , 



(6) 

0.25 GeV-2. 
= 13.4, and 



/j/ij, = 11.2. The other parameters in Eq.(j4]) have been determined by fitting[27[ the total 
cross section data of the photo-production of p and u : Pu = Pd = 2.07 GeV~^ and /Iq = 1.1 
GeV^. 

With the parameters specified above, our task is to examine the extent to which the 
total cross section of photo-production of J/\l/ can be fitted by only adjusting the Pomeron- 
charmed quark coupling constant Pc- This will be discussed in section IV. 
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B. Pion-exchange amplitude 

We observe from Particle Data 26| that the width of the J/^ — Tr^p^ is significant, 

Tj/^^^opo = 0.92 MeV x (0.56 ± 0.07)% . (7) 

With the vector meson dominance (VDM) assumption, this experimental information allows 
us to calculate the one-pion-exchange amplitude of 7 + — t- J/\& + iV, as illustrated in 
Fig{T]|^b), by using the following Lagrangian 

L = Lj/\j,^p0^o + Lt^nn + LvDM (8) 

with 

LnNN = -^^^^Nlf.lbT^Nd'' ■ (pn, (10) 

LvDM = ^^A^pI , (11) 
Jp 

where p°, (f)j/';j,i3, (f-K, and i/jn are the field operators for p°, J/\l/, vr, photon (7), and 
nucleon (A^), respectively. The mass for particle a is denoted as m^. The well determined 
coupling constants are f^Niy/4:7i = 0.079, e^/47r = 1/137, and fp = 5.33. To determine 
gj/'b,pOTTO, we use Lj/^^pO^o given in Eq.(|9]) to calculate the decay width 



{2n)\ V / dQ,\ < kXp\H\p = 0, Xj/^ > \2 kEAk)Ep{k) 



where k is defined by rrij/^ = ET,{k) + Ep{k), and 

1 1 



X 6^'^"^fc^6.,,^ (fc^P.^/^.A,,, (P) [ ^ -.^ "^"f^ ^ ]^ . (13) 



Here we have included a dipole cutoff function with a range parameter Aj/^. The four- 
momenta are defined in the rest frame of J/\l/: 

p = (mj/q,,0) , 
F = (i?p(A;),^), 
r = (K(A;),-^). 

By using Eqs.f fT2|) -f fT3|) and the experimental value given in Eq.(I7]), we find gj/qi^nOpO = 0.032 
for a cutoff Aj/^ = 2000 MeV. 

With the Lagrangian Eq.®, the one-pion-exchange invariant amplitude for 'y{q)+N{p) — >• 
J /'^{k) + N{p') can be written as 

If^ = Un^dP'Kik, Aj/^)Mr(p',p)e.(g, X^)u^M (14) 
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with 



Mr(p',p) 



G X F{t)-^e^''^^Kqp[j ■ ip' - p)] , 
t — mi 



(15) 



where t = {p — p'Y, and 



F{t) = FT,NN{'t)Fj/q,^p0^o{t) . 



(16) 



(17) 



Here we have introduced a cutoff form factor F{t) to regularize the interaction verteces. For 
simphcity, we use the following form 



We set n = 4, and A = Aj/^ = 2000 MeV. 

III. PHOTO-PRODUCTION OF [B]j/^ BOUND STATE 
A. Reaction Mechanism 

With the impulse approximation, we assume that a J/\l/ is produced on a nucleon in 
the target nucleus A and then is attracted by a spectator system B to form a bound state 
[Sjj/,^. For simplicity, is denoted as d in the following formula. 

With the mechanism illustrated in Figj2l the cross section of '~f{(l)+A{—q) — > N{p)+d{—p) 
in the center of mass system can be written as 




(18) 



da 

dn 



{2TrYEA{q) pEN{p)Ed{p) 

w w 




(19) 



Y 




N 



A 



FIG. 2: The impulse approximation mechanism of 7 + A — )■ + [S] j/if reaction. j4 is a nucleus 
with mass number A and B could be a nucleus with mass number {A — 1) or a [g^^"^"^)] multi-quark 
cluster. 
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where 



X 



(20) 

Here a]„_^^ is the creation operator for a J/"^ with wavefunction Xj,mj{Qd,mj/xi,), bj^^mj^ an 

annihilation operator for a nucleon with wavefunction (pj^'maiQAi^sa)-, and the 7 + — >■ 
J/'^ + N amphtude is 



^ - / U/TTAI-\ - ^ - - , rriN i TUN 

< puis, kmj/q,\t{W)\qX,Pa'ms^ > 



(27r)3 ^2Ej/^{k)\j ENiPa) v^V E^ip) 

x[UraAPK{k,mj/^){Mi:''{p,p^) + M;''{p,Pa)} 

(21) 

where M^" is the pion-exchange amphtude given in Eq.f llSp . and A^p'' is the Pomeron- 
exchange amphtude in Eq.Q. 

For simphcity, we wiU only perform calculations for the reactions on ^He and ^He. For 
estimations of cross sections on these target nuclei, it is sufficient to use the s-wave harmonic 
oscillator wavefunctions for both the target A and B in the [-B]j/>ir bound state. We also 
only consider the case that the J/\E' in the produced bound -Bj/* is on an s-wave orbital. 
For the case B = A — 1 nuclear system, we thus write the initial (1$"^ >) and final (|\E' >) 
nuclear states as 

|$^> = [\N > ^\<!>^-^ >]l=o , (22) 
1^ > = > >]^^o, (23) 

where L is the relative angular momentum between or J/\& and the {A — 1) nucleus. 
Explicitly, we have 

I'^J^M.j, > = J2 JdMd\jJA-im,Mj^_^ > sLj'^'j;i,Af,^_^ > • (24) 

Then the momentum variables in Eqs. fl20p and fl^T]) are 

Pa = p + k-q, (25) 
P/3 = -p - k = -q- Pa, (26) 

^ _ PaEA-ljpp) -ppEN{Pa) 

Ea-,{P,) + E^iPa) ' ^''^ 

Qa = -PP^J/"^^^) (28) 

Ej/^(k) + Ea-i{pp) 



where Qd {Qa) is the relativistic relative momentum between J/^ (A^) and the {A — 1) 
nuclear system. 

For the target ^iJe, we have J^i = and assume that m > the '^He ground 

state with Ja-i = 1/2. We then have the following simplicities: 

\^^A,Mj^ > = < JAMj^\]^JA^imj^Mj^_^ >< JdMdljJA-irrijMj^^^ > 

^ ^ < JdMjJjaJ - nij^mj > , (29) 



and 

1 



— * 

-T^nQd) , (30) 



</'i.,rn,.J<5A,msJ = 5j,,i/25mj„,m,^^=i?(gA) , (31) 

v47r 



Eq.( l20|) then becomes 



M, > 

= Yl Yl < JAMjJjaJA-imj^Mj^_^>< JdMd\jJA-imjMj^_^> 

/■ - 1 - 1 

x[ / dk^=F{Qd) < pmskmj/^\t{W)\qX,Pams^ > ^=R{Qa)] ■ (32) 
J v47r V47r 



We have applied the formula Eqs.( IT9|) and (!32|) to estimate the production cross section 
on ^ife. We use the usual s-wave harmonic oscillator wavefunction with b = 1.32 fm for the 
target "^He 

R{p) = [Ne^] (33) 

with the normalization J R^{p)p'^dp = 1. For the J/"^ wavefunction in c? = pife]j/^, 
we will generate a s-wave ipj/'S'i''^) from a potential Vj/t^sif') with the normalization 
J r'^dr\iljj/ii,{r)\'^ = 1. The wavefunction in Eq. fl30|) and also in (1321) can then be calculated 
from 

POO 

F{p) = / r'^drjo{pr)iljj/<!,{r) , (34) 
Jo 

where jo{z) is the spherical Bessel function. The form of Vj/,i,^B(r) will be discussed in 
section IV. 

The above formula can be easily extended to investigate other possible impulse approxi- 
mation mechanisms as far as all wavefucnctions in the bound A and [5] j/^ are all in s waves. 
This is what we will need in section IV when we consider the production of [q'^jj/* from the 
g^-N component of ^He. 
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IV. RESULTS 



A. Models of 7 + iV ^ J/'I' + reaction 



We first develop a model consisting of Pomeron-exchange and pion-exchange mechanisms, 
as described in section II. With the parameters specified there, we try to fit the available 
total cross section data of 7 + p — )■ J/\l/ + p up to invariant mass W = 300 GeV by only 
adjusting the charmed quark-Pomeron coupling constant Pc- With /3c = 1.21 we only able to 
fit the data up to 20 GeV. Clearly, the result at high energy is not satisfactory as shown in 
the red dashed curve in the left-hand side of FigjSl We then find that by changing ao of the 
Regge trajectory in the Pomeron propagator Eq.([6]) from ao = 1-08, as determined in the 
previous fits|27| to the total cross sections of p and u photo-production, to ao = 1-25, we 
are able to get a very good fit to the data by choosing (3^ = 0.84 GeV~^. Our fit is the solid 
black curve in the left-hand side of FigJSi We thus will use the model with ao = 1.25 and 
/3c = 0.84 GeV~^ (PM model) in our investigations. As also seen in the insert in the left-hand 
side of FigJ31 the contribution (magenta dotted curve) from the pion-exchange amplitude, 
as defined by Eqs. (fT4|) - (fT8|) . is very weak except in the very near threshold region. 
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FIG. 3: The total cross section of the 7 + p — )• J/V' + P reaction as function of the 7p invariant 
mass W. Left: the black solid and red dashed are the results from using (ao = 1.25, Pc = 0.84 
GeV~^; called the PM model) and (ao = 1-08, Pc = 1-21 GeV~^) within the model which include 
both the TT-exchange and Pomeron-exchange mechanisms. The magenta dotted curve in the insert 
is the contribution from the vr exchange. Right: the black solid, blue dotted and green dot-dashed 
lines are from the PM model, th e 2q model of Ref. 13], and the 2g + 3g model based on Eq.(| 
The experimental data are from 291437 [. 



We next consider the model of Ref. l2] based on the two-gluon {2g) and three-gluon {3g) 
exchange mechanisms. In terms of the normalization defined by Eq.(l2]), the amplitude of 
this model can be written as 



(27r)3 ^2Ej/4k)y E^{k) V E^{q) 
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(35) 



with 



= ^^e^*/^ (36) 
M3. = ^^^e^*/^ (37) 



2 



2mNmj/^ + mj,< 

X = , . 38 

- mi 



where R = 1 fm, h = 1.13 GeV ^ are taken from Ref. [1^. We follow Ref. 12| to determine 
the parameters and A^g by fitting the data up to only 20 GeV. In the two-gluon-exchange 
model {2g), we set = and obtain = 0.028 MeV~^ from the fit. In the 2g + 3g 
model, the fit is obtained by choosing A2g = 0.023 MeV"^ and Asg = 2000 MeV'^ The fits 
for the 2g and 2g + 3g models are the dotted and dot-dashed curves in the right-hand side of 
FigJSl respectively. Clearly, they have differences with that (black solid) of the PM model, 
as can be seen more clearly in the insert in the right-hand side of Figj3l Here we also see that 
the data in the region near the J/\E' production threshold are very limited and uncertain. 
We will therefore perform calculations using the PM, 2g, and 2g + 3g models to examine 
the model dependence of our predictions. Clearly, precise data in the near threshold region 
are needed to make progress. 



B. Photo-production of J/^-Nucleus bound states 

Following the previous investigations 0, [28[, we assume that the interaction between a 
J/\E' and a nucleus with mass number A can be parameterized as a non-relativistic potential 
of the following Yukawa form 

K7/*,A(r) = -ttA^— . (39) 

There exists two different approaches to determine the parameters aA and fiA for the nucleon 
with A = 1. We will explain these in the following two subsections. 



1. Pomeron- quark coupling model 

Motivated by the previous studies in Quantum Electrodynamics, it is assumed in the 
approach of Ref.[3] that the J/\i/-A forward angle scattering amplitude at very high energy 
can be related to the matrix element of the potential Eq. (!39|) which is understood to be 
valid only in the region where J/\i/ moves non-relativistically. They further assume that the 
J/\i/-A amplitudes can be calculated by using the Pomeron-exchange model of Dannachie 
and Landshoff jsj. In the very high energy approximation, the differential cross section of 
J/\E'-A elastic scattering can be related to the parameters a a and fiA of the potential Eg. (139 p 
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TABLE I: Parameters for determining the potential Eq. (j39|) using the Pomeron-quark couphng 
model defined by Eqs. (j4ip - (j43p . The predicted binding energies (B.E.) for proton (A = 1), ^He 
(A=3) and (A=12) are also listed. 



A 
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P 
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(GeV 




(iVleV j 


1 


3.9 


0.63 


1.85 


1.85 


0.64 










2.05 


1.21 


0.47 










2.05 


0.84 


0.33 




3 


9.5 


0.26 


1.85 


1.85 


0.33 


19.86 








2.05 


1.21 


0.23 


3.27 








2.05 


0.84 


0.16 


0.04 


12 


12.69 


0.19 


1.85 


1.85 


0.73 


280.0 








2.05 


1.21 


0.53 


165.0 








2.05 


0.84 


0.37 


67.0 



by following relation 



|(7/v. A J/^ A) ^ P^.^^/.(^)]-[3^/^V^^^(^)]- (40) 

- (41) 



where t is the momentum-transfer squared, Pu/d if^c) is the Pomeron coupling with the 
up/down ( charmed ) quarks, Fji^,{t) and FA{t) are the form factors for J/\i/ and the 
nucleus with mass number A, respectively. They further assume that in the t — ?■ limit, the 
slope of da/dt is mainly determined by dFA(t)/dt and that FA{t) can be identified with the 
nuclear electromagnetic form factor. One then gets the following relations 

f^A - l^^l*=o-— (42) 
aA = ^ f^A- (43) 



The radius < R\ >^/^ can be taken from Ref. 38|. The Pomeron-quark coupling constants 



can be taken from fits to the data of meson-nucleon scattering or photo-production of vector 
mesons. Once a a and fj,A of the potential Eq.( l39|) are determined, we can predict the possible 
J/\l/-nucleus bound states. In Table [H we list our results for proton {A = 1), ^He (A=3) 
and ^^C (A=12) for various sets of Pomeron-quark coupling constants. The first rows in 
the results for each A are based on the flavor independent (S^/d = Pc = 1-85 GeV~^ of 
Ref.0]. The other two results use the Pomeron-quark coupling constants Pu/d = 2.05 GeV~^ 



determined [27| in the fits to the data of photo-production of p and cu, and Pc determined 
from the fits described in section IV. A. 

With the determined potential parameters and ^a, the predicted binding energies 
(B.E.) for each considered nuclear system are listed in the last column of Table [H For the 
A = 1 case, we see that there is no J /'^-N bound state. But all three models predict bound 
^He]j/^, and states. In the left-hand side of FigJU we show the predicted cross 
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sections of 7 +^ He — )■ pifejj/^ + n. Clearly, the predicted cross sections depend on the 
Pomer on- quark coupling constants. Furthermore, their magnitudes depend sensitively on 
the binding energy (B.E.) of the predicted [^ife] j/^ system. As the binding energy decreases 
from 19.86 MeV to 0.04 MeV, the predicted cross sections drop by two orders in magnitude. 
This can be understood from the right-hand side of FigJH where we compare the J/^-^He 
relative wavefunctions which are used in predicting the cross sections in the left-hand side. 
We see that the wavefunction (solid black) for B.E. = 19.86 MeV is much shorter range than 
the other two cases and hence gives more cross sections in this large momentum-transfer 
reaction. This is explicitly illustrated in Figj5] where we show that the cross section (red 
dashed curve) calculated from keeping only the high momentum part (pj/^ > 1400 MeV) of 
the J/'^ wavefnction in the integration in Eq. fl32|) is very close to the full calculation (solid 
black curve). 

In FigEl we see that the predicted differential cross sections are forward peaked, as 
expected from the Pomeron-exchange mechanism. In FigjTl we show that the predicted cross 
sections depend on the 7-l-A^ — )■ J /'^ + N model. Their maximum values are, however, 
comparable ~ 0.1 — 0.3 pico-barn. Clearly, it is important to get accurate data of 7 + iV — )■ 
J /'^ + N aX low energies to refine the employed model for making more precise predictions. 

The P^C] j/^, can be produced by 7 C [^^C] J/* + f^- However, making predictions 
for the cross sections of this process is beyond the scope of this paper since the simple s-wave 
description of the nuclei in section III is no longer a reasonable approximation for nuclei 
heavier than '^He. 




W(GeV) R(fm) 



FIG. 4: The total cross section of 7 + \^He] — t-^ He^j/^^ +7i as function of 7-^-ffe invariant mass W 
(left) and the wave function (right) for J/ip — ^ He system. The black solid, red dashed and blue 
dotted-dashed lines are calculated by using the potential Eq. ([39|) with A = 3, jiA = 0.257 GeV and 
aA = 0.33,0.23,0.16, respectively. 



2. Folding model 

While all three J/'^-N models listed in Table [T] do not have bound states, there exist a 
possibility that adding the J /'^-N interactions from the nucleons in a nucleus could lead to 
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0.4 




W(GeV) 

FIG. 5: The cross section of 7 + [^He] — )• ['^i?e][j/^] + re as function of invariant mass W. 

The red dashed curve is obtained from keeping only the contribution from the J/^ wavefunction 
with k > 1400 MeV in the integration of Eq.([32]). 




30 60 90 120 150 180 
e (Degree) 



FIG. 6: The differential cross section of 7 + [^He] — t- pi?e][j/^] + re vs the angle of out going N 
with the center mass 6.9 GeV. 

bound states. To explore this possibility, we follow the usual nuclear physics approach to 
construct a folding potential for the interaction between a J/\l/ and a nuclear system, 

K//*,A(r) = j vj/q,^N{r- f')pA{r')df' , (44) 

where wj/^,Ar(r) = Vj/^_i(r) as defined by Eq. (!39l) with A = 1, and the nuclear density is 
normalized by 

j pA{r')dr' = A. (45) 

For ^He we use pA(r) = po^"'^^^^ with h = 1.32 fm which is obtained by fitting the ^He 
charge form factor at low momentum-transfer. For heavy nuclei, we use the Woods-Saxon 
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FIG. 7: The total cross section of 7 + [^i^e] — )• pi?e][j/^] + n vs the certam mass of system. The 
black solid, red dashed and blue dot-dashed lines are for the pomeron and vr exchange, Brodsky's 
2g model and 2g+3g model, respectively. 



form 3^ 



p^^"^ = p \ + el-^y^ ^^^^ 



with R = g^nd t = 0.53 fm. 

Our results using the parameters of Ref. [3] to calculate vj/<i,^M{R) in Eg. (HI]) are listed 
in the first row of Table [TTl We see that the folding model gives 1.62 MeV (7.0 MeV) for 
[^He]j/^ ( [^^C]j/^) which are much less than 19.86 MeV ( 280 MeV) listed in Table E The 
predicted cross sections for — )■ ^He\j/^, +n are also found to be much weaker, close 
to the blue dot-dashed curve (0^3 = 0.16) in FigJH Clearly, it is difficult to measure such a 
loosely bound piJe]j/^ state. 

To examine the model dependence, we also consider folding potentials by using three 
other J/'^-N models. Two0,l6[ of them are constructed by using the results from the heavy 
quark effective field theory calculation by Peskin[l|. The third onejsf is from Lattice QCD 
calculation. Their results can also be written in the Yukawa form of Eq.f lS^ with A = 1. We 
find that these three models do not generate a [^He]j/ifi bound state as indicated in Table 
im For ^^C, the binding energies from folding model are much weaker than those listed in 
Table |T] from the Pomeron-quark coupling model. 



C. Photo-production J/^-{q^) bound states 

In Ref . [ill , it was suggested that a cc system could interact strongly with the color octet 
3-quark [g^jg component of the six-quark cluster ([g^] = [g^jsf^^Js) which could dominant the 
short-range part of the deuteron wavefunction. The possible attractive force between a J/^ 
and a six-quark cluster was suggested in the study of Ref. (Tlj where the excitation of a hidden 
charm \qqqqqqcc > state is introduced to explain the spin correlation oi pp elastic scattering 
near the J/'^/ production threshold. Here we examine the condition under which a bound 
[q^]j/^ color singlet state can be produced in the 7 +^ He — )■ [g^]j/ij + N reaction. Unlike 
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TABLE II: The binding energy of J/?/'- nucleus calculated with the folding potential defined by 
Eq. (|44p with parameters of Vj/x^ j\f = —ai^—y— taken from different references. The parameters 
of Ref.Q] are obtained from reproducing the scattering length a = —0.24 fm given in Ref.Q]. (A 
Gaussian form of the Vjj^,^^ was used in Ref.0]) 

Model Parameter (MeV) Binding Energy (MeV) 
ai ^1 (GeV) [H]j,^ ?He]j,^ PC]j,^ 

Ref.[a 0.64 0.63 - 1.62 7.0 

Ref.[4j 0.20 0.63 - - 0.91 

Rcf.[5] 0.10 0.63 - - 0.003 

Rcf.[6] 0.06 0.63 



the predictions for the photo-production of pife]j/^ described in the previous subsection, 
very little information on [q^] and the [q^] — J/^ interaction is available. We thus need to 
make various assumptions which can only be considered to be plausible for estimating the 
production cross sections. 

In the impulse approximation, as described in section II, we need the initial A^-[g^] 
wavefunction in ^He and the final J/'^-[q^] wavefunction to calculate the cross section of 
7 +3 He — [q'^] J/* + N . In the following subsections, we explain our procedure for modeling 
these two ingredients of our predictions. 



1. Wavefunction of N -q^ 



6 3 



He 



We start with a formulation of Refs. 17, 41 within which the Hamiltonian for a two 



nucleon system is written as 



H = Ho + Vjyjy + ^ h 



(47) 



where a denotes collectively the total angular momentum J, the total isospin T, and the 
parity P, and v^n is a meson-exchange nucleon- nucleon (NN) interaction. The vertex 
interaction h^gSja^^^ defines the formation of a six-quark state [g^]" in A^A^ collisions. The 
six-quark states [q^]" are identified with the states predicted by the Bag model calculations 
of Mulder[l3]. By appropriately choosing the form of the vertex interaction /i[g6]a^7VAr, 
the NN scattering amplitudes derived from the Hamiltonian Eq.( l47|) are identical to those 
given by using the P-matrix approach of Jaffe and Low 13|] and the Compound Bag Model 



formulation developed in Refs. |l5l . Il6 



We will make use of the results of Fasano and Lee|l7l. Il8|. They determined the mass 



of [g^]" cluster and the interaction h[q6^c.^^^ for a =^ Sq and "^Si by fitting the A" A" scattering 
phase shifts up to 1 GeV. Within the simple s-wave harmonic oscillator model for ^He, the 



probabilities P[qG]a of finding the [q ]°'-N in He are estimated[l8| to be P, 



and P[q6fsi 

of [q^] determined in Ref. 



[g6]iSo 



0.7% 



0.06%. For simplicity, we neglect the small ^Si component. The bare mass 

TtI is Miso = 2150 MeV. Here we will use these information to 
model the relative wavefunction of [q^]'^-N which is needed to calculate the cross sections of 
7+3 He^[q%/^ + N. 
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R(fm) 

FIG. 8: The solid curve is the normahzed density distribution calculated from using Eqs. (j50p - (|5ip 
with b = 1.35 fm and Ic = 0.5 fm. The dashed line defines Tc = 0.29 fm for obtaining the probability 
Pr gnigg = 0.7% for finding the [q^] — N component in ^He by using Eq.([53|). 



We assume that the charge distribution in the region with the distance r < from 
the center of ^He is completely due to [q^Y'-N components of the wavefunction. This is 
illustrated in FigJHl Each Tc clearly corresponds to a choice of P^e. Within such a model, 
the charge form factor of ^He is written as 

Fcm = F^\Q') + Ff-''{Q% (48) 



We next observe that within the conventional nuclear model[40|, the impulse approximation 
(lA) calculation, which includes only the one-body nucleon current, of -Fc(Q^) is very close 
to the data in < about 10 fm~^ and can be reproduced very well by the Gaussian 
distribution of the s-wave harmonic oscillator wavefunction. The lA results from Ref. 





in this region are the solid squares in FigJH We next demand that the s-wave three- 
nucleon wavefunction reproduce these lA results. In addition, the resulting Fc(Q^) in the 
higher region must have the similar structure of lA up ~ 20 fm~^ although we do not 
have higher partial wave components of the three-nucleon wavefunction. We achieve this by 
using the s-wave harmonic oscillator wavefunction with Jastrow two-body correlation used 
in Refs. 0,1431. We write 



Fc{.Q') = I e~'^-'p{r)dr (49) 

with 

p(fi) = / df2p2iri,f2), (50) 



where the two-body density is defined by 

p2{ri,f2) = Ne 26^^(1 



2 I 2 
''l+'-2 
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^ (V^6)3\/ [Z2 + 62]3/2_/3- (51) 

As seen in FiglH the solid black curve calculated with b = 1.35 fm and Ic = 0.5 fm can 
reproduce the impulse approximation calculation results (solid squares) given in Ref. 401] up 
to ~ 10 fm~^. At higher Q^, the solid curves have the similar structure of the lA results. 
For our present s-wave calculations, we consider the solid curves in Figl9]as the Fc{Q'^) in 
Eq. lHHl) . Accordingly, the NNN contribution in Eq. fl48l) is calculated from 



/oo 
r^dr I rfl],e"'«->(f) (52) 

and the probability P[g6] is defined by 

P^6 = / r'dr I dVtrpir) . (53) 



For PqH = 0.7% determined in Ref.[l7| within the Compound Bag Model of NN scattering, 
we choose Tc = 0.292 fm to calculate Eq.f l52p and get the blue dotted curve in the left-hand 
side of Figjni In the right-hand side, the blue dotted curve is from the calculation using 
Eq.f l52|) with Tc = 0.630 fm which gives P^e = 6.3%. Clearly, both results agrees well with 
the lA (solid squares) and the solid curve only in the low region. Our next task is to 
model F'^ ~^(<5^) such that for each Vc, FdQ"^) (solid black curve) in Figl9]up to ~ 15 
fm^ can be reproduced, as required by Eq. (l48!) . 




FIG. 9: The charge form factor FJ Q"^) for ^He. The solid squares are from the Impulse Approx- 
imation (lA) calculation of Ref. 40|. The dotted blue curves are F^^ {Q"^) calculated from using 



Eq.([52]) with rc = 0.292 fm and P^e = 0.7% (left) and = 0.630 fm and P^e = 6.3% (right). The 
red dashed curves are from adding the — N contributions calculated from using Eq. ()54p with 
h* = 0.185 fm (left) and b* = 0.414 fm (right). 



For simplicity, we assume that Ff-^{Q^) can be calculated from a normalized Gaussian 



distribution 



Ff-^{Q) = P,6 / e~'^^'[y4^ew]df . (54) 
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Accordingly, the mean radius of q^-N can be defined by 



^2 



< >= / h r^u^Yi ^'^Vdr. (55) 



We adjust h* for each Pqd to fit the sohd curves in FigJHi We find that if P^e is larger 
than 6.5%, no h* can fit the form factor defined by Eq.f l48l) in the < 20 fm~^ region. In 



the cluster model of Ref.[20||, P„6 = 4.0% is obtained from fitting the He form factor. The 
Pg6 = 15% determined in Ref. 19| by fitting the structure function of ^He{e, e') is beyond 



what our formulation can accommodate. For comparison, we thus choose three different 
models with P^e = 0.7%, 4%, and 6.38% for our calculations. In Table UTTl we list Vc, b*, 
and also < >^/^ calculated from using Eq. (l55|) for these three cases. Our fits are the red 
dashed curves in Figl9]for P^e = 0.7% (left) and 6.38% (right). 

Once b* is determined, we then assume that the relative wavefunction of [g^]-N can be 
described by the harmonic wavefunction with the same b*. This should be reasonable for 
making order of magnitude estimates in this work. A more sophisticated approach should 
account for the quark charge distribution in [q^] which is beyond the scope of this work. 
Also, the sharp cutoff at r = Tc to define F^^ in Eq. fl52l) should perhaps be better modeled. 
For our present qualitative estimations, this simple procedure should be sufficient. 



TABLE III: The parameters for ^q-N and ^q-J/tp systems. See text for the explanations of the 
notations. 



q^-N 




^-J/iP 


Model 


Pg^ Tc b* 


< >V2 


flqfi 


OgG B.E. 




(fm) (fm) 


(fm) 


(GeV) 


(MeV) 




0.7% 0.292 0.185 


0.226 


0.6 


1.33 498.42 








1.0 


1.50 389.63 


Bi 


4.0% 0.533 0.346 


0.424 


0.6 


0.83 104.08 


B2 






1.0 


1.05 79.77 


Ci 


6.38% 0.630 0.414 


0.507 


0.6 


0.75 65.84 


C2 






1.0 


0.97 51.33 



2. Wavefunction of ^q- J ftp bound state 

We follow the procedure of subsection IV. B to assume that the q^-J/^ bound states 
([g^] j/ij,) are also defined by a potential, of Yukawa form 

yj/^i-.q^ir) = -Oiqs— — . (56) 

We expect that if a bound state can be produced, its size must be small for color 

field to give strong attractive force. Thus it is reasonable to assume that the mean radius 
of [g'^jj/-!- is close to the value < >^/^ of the initial q^-N system listed in Table UTTl We 
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find that such a small size can be generated from choosing /i^e > 0.6 GeV in defining the 
potential Eq.( l56|) . Once a value of /x^e is chosen, we then determine the potential strength 
ag6 by requiring 

<r'>= I |0,a,,/^(r)|Vrfr, (57) 

where 0q6 j/,j,(r) is the J/'^-q^ relative wavefunction generated from the potential Eq.f l56p . 
and the values of < >^/^ for various considered cases are listed in Table lllli The re- 
sulting ag6 and the binding energies (B.E.) are also listed there. Here we note that the 
binding energy increases as the mean radius < >^/^ and the corresponding probability 
Pqd decrease. 



3. The Results of -y+^He^ [q^].j/i, + N 

With the wavefunctions for q^-N and q^-J /'^ specified in the previous subsections, we can 
use the formula in section III, with trivial changes in notations and spin quantum numbers, 
to calculate the total cross section of 7 +^ He [q^]j/i> + N . However, we need to multiply 
the results by the probability PqH of the A^-[g^] component in ^He\ namely the results from 
using Eq. (fT9l) is changed to 

d(J „ fin ,da ^ 

where [^]o is calculated from using Eq. ffTOj) and all subsequent equations in section HI. A. 

We first consider the case that Pj^s] = 0.7% as determined in Refs. 17|, [iSj from fitting 
the A^A^ phase shifts up to 1 GeV. By using the parameters for models Ai and A2 listed 
in Table IIIIl we obtain the results shown in FigJTOl We observe that with the same small 
radius < >i/2= 0.226 fm for the produced system, the predicted cross sections are 

very close despite their potential range, measured by l//ig6, and couphng constant a^e can 
be very different. The same finding is also from comparing the predicted cross sections from 
the models Bi and i?2, and also the models Ci and C2. 

In the left-hand side of FigJUl we show the dependence of the predicted cross sections 
on Pqd by comparing the cross sections from three models Ai, Bi, and Ci listed in Table 
mil We observe that as P^e decreases, the peak is shifted to higher energies. Each case 
has different threshold energy due to their differences in binding energies, as seen in Table 
mil Their magnitudes are comparable despite their P^e are very different. We find that 
this is due to the fact that the cross section [^]o in Eq.(l58l) for the model with smaller 
Pq6 = 0.7% is a factor of about 10 larger than that for the model with larger P^e = 6.38%, 
since this large momentum transfer reaction favors the production of [q^]j/m with smaller 
size characterized by < >^/^ in Table lllli The situation is similar to what we discussed 
in explaining the results shown in FigJH Thus the magnitudes of the cross sections from 
three models at peak positions are comparable because the factor of about 10 difference in 
[^]o in Eq.( l58l) is compensated by the similar factor of about 10 in P^e. However, the three 
models have rather different energy dependence, as also seen in the left-hand side of FiglTTl 
On the other hand, they are all forward peaked, as shown in the right-hand side of FigJTTl 
for the differential cross sections aXW = 6.6 GeV. 

The results shown in Fig{TT] suggest that the upper bound of the predicted total cross 
sections of 7 -1-'^ He — )■ [q'®] j/<t + A^ is about 2-4 pico-barn 
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FIG. 10: The total cross section oi ^ +^ He — )• [q^]j/ti, + N as function of the j-^He invariant mass 
W. The black solid and red dashed curves are for the case Ai and A2 in the Table HIH respectively. 
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FIG. 11: The black solid, red dashed and blue dotted-dashed lines are from Models Ai, Bi, and 
Ci of Table HIH respectively. Left: The total cross section of 7 +^ He — )■ [q^]j/^ + iV as function of 
the 'y-^He invariant mass W. Right:The differential cross section of 7 +^ He — )• Q[j/^] + of the 
outgoing N at ^-^He invariant mass W = 6.6 GeV. 



V. SUMMARY AND DISCUSSIONS 



We have presented predictions of the cross sections of 7+^ife — >■ iV+ pife] j/* reaction at 
energies near the J/\l/ production threshold. In the impulse approximation, the calculations 
have been performed by using several 7 + A^— > J /'^ + N models based on the Pomeron- 
exchange and pion-exchange mechanisms. The J/^ wavefunctions in pife] j/^ are generated 
from various J/\l/-nucleus potentials which are constructed by either using a procedure based 
on the Pomeron-quark coupling mechanism or folding a J/\E^-N potential Vj/x^i^n into the 
nuclear densities. We consider vj/ifi^N derived from the effective field theory approach. Lattice 
QCD, and Pomeron-quark coupling model. The upper bound of the predicted total cross 
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sections is about 0.1 - 0.3 pico-barn. 

Clearly, our investigations are only for estimating the cross sections to facilitate the 
experimental considerations of possible measurements of piiTejj/^ bound states at Jefferson 
Laboratory. Several improvements are needed for more quantitative predictions. First we 
need precise data of7 + p^' J /'^ + p near threshold to distinguish several models we have 
considered and also to develop a more sophisticated model. We also need the data to pin down 
the J/\Ef-N interaction for a more realistic calculation of J/\Ef-nucleus potential such as the 
folding model considered in this work. One possibility is to use the vr+^if— )-p + J/\& + n 
reaction to extract the J /'^-N scattering length, as suggested in Ref.jl]. Alternatively, 
we can apply the model presented in this paper to determine the J/'^-N interactions by 
investigating the — p+J/'^+n reaction. Possible experiments on these two processes 

will be very useful. We of course also need to use more realistic wavefunctions for ^He and 
^He while the s-wave oscillator wavefunctions employed in this investigation are reasonably 
consistent with the charge form factors calculated from the conventional nuclear models. 

Motivated by the previous investigations ll|, |l2| on the effects due to multi-quark clusters 
in pp and 7 -|-^ if — )■ J /'^ + n + p, we have also considered the possibility of the production 
of a [q^] — J/^ bound state due to a six-quark [q^] cluster in ^He. The Compound Bag 
Model of A^A^ scattering and the quark cluster model of nuclei are used to estimate the 
[g^]-N wavefunction in ^He by imposing the condition that the sum of the contributions 
from [g^]-N and NNN components to the ^He charge form factor must be consistent with 
what are predicted by the conventional nuclear models 40| which explain the data very well. 
The upper bound of the predicted total cross sections of 7 +^ He — )■ + A^ is about 

2-4 pico-barn, depending on the model oi ''y + N J /'^ + N used in the calculations. If 
such bound states can be identified, it will open up a new window for investigating the role 
of the gluon field in determining the hadron structure. 
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